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Abstract —This paper considers a scenario of short-range com- 
mnnication, known as device-to-device (D2D) communication, 
where D2D users reuse the downlink resources of a cellular 
network to transmit directly to their corresponding receivers. 
In addition, multiple antennas at the hase station (BS) are used 
in order to simultaneously support multiple cellular users using 
multiuser or massive MIMO. The network model considers a 
fixed number of cellular users and that D2D users are distributed 
according to a homogeneous Poisson point process (PPP). Two 
metrics are studied, namely, average sum rate (ASR) and energy 
efficiency (EE). We derive tractable expressions and study the 
tradeoffs between the ASR and EE as functions of the number 
of BS antennas and density of D2D users for a given coverage 
area. 


I. Introduction 

The research on future mobile broadband networks, re¬ 
ferred to as the fifth generation (5G), has started in the past 
few years. In particular, stringent key performance indicators 
(KPIs) and tight requirements have been introduced in order 
to handle higher mobile data volumes, reduce latency, increase 
connectivity and at the same time increase energy efficiency 
(EE) 0, 0. The current network and infrastructure cannot 
cope with 5G requirements—fundamental changes are needed 
to handle future non-homogeneous networks as well as new 
trends in user behavior such as high quality video streaming 
and future applications like augmented reality. 5G technology 
is supposed to evolve existing networks and at the same time 
integrate new dedicated solutions to meet the KPIs 0. The 
new key concepts for 5G include massive MIMO (multiple- 
input multiple-output), ultra dense networks (UDN), device-to- 
device (D2D) communications, and huge number of connected 
devices, known as machine-type communications (MTC). The 
potential gains and properties of these different solutions have 
been studied individually, but the practical gains when they 
coexist and share network resources are not very clear so far. 
In this paper, we study the coexistence of two main concepts, 
namely massive MIMO and D2D communication. 

Massive MIMO is a type of multiuser MIMO (MU-MIMO) 
technology where the base station (BS) uses an array with 
hundreds of active antennas to serve tens of users on the 
same time/frequency resources by coherent precoding 0.0- 
Massive MIMO techniques are particularly known to be very 
spectral efficient, in the sense of delivering sum rates 0. 
This comes at the price of more hardware, but the solution 


is still likely to be energy efficient & 0 - On the other 
hand, in a D2D communication, user devices can communicate 
directly with each other and the user plane data is not sent 
through the BS 0. D2D users either share the resources 
with cellular networks (overlay approach) or reuse the same 
spectrum (underlay approach). D2D communication is used for 
close proximity applications and can bring tremendous gains 
in terms of offloading core networks and achieving higher data 
rates with less transmission energy. 

We consider two network performance metrics in this work: 
The average sum rate (ASR) in bit/s and the EE defined as the 
number of bits transmitted per Joule of energy consumed by 
the transmitted signals and the transceiver hardware. It is well 
known that these metrics depend on the network infrastructure, 
radio interface, and underlying system assumptions 
Goi. The motivation behind our work is to study how the 
additional degrees of freedom resulting from high number of 
antennas in the BS can affect the ASR and EE of a multi-tier 
network where a D2D tier is bypassing the BS. We focus on 
the downlink since majority of the payload data and network 
energy consumption are coupled to the downlink 0. We 
assume that each D2D pair is transmitting simultaneously with 
the BS in an underlay fashion. In addition, we assume that the 
communication mode of each user (i.e., D2D or cellular mode) 
has already been decided by higher layers. 


A. Related Work 


The relation between the number of BS antennas, ASR and 
EE in cellular networks has been studied in 0, 0, IID^ 

1121 among others. The tradeoff between ASR and EE was 
described in 0 for massive MIMO systems with negligible 
circuit power consumption. This work was continued in Ol 
where radiated power and circuit power were considered. In 
0^ joint downlink and uplink design of a network is studied 
in order to maximize EE for a given coverage area. The 
maximal EE was achieved by having a hundred BS antennas 
and serving tens of users in parallel, which matches well 
with the massive MIMO concept. Eurthermore, the study | [T^ 
considered a downlink scenario in which a cellular network 
has been overlaid by small cells. It was shown that by 
increasing the number of BS antennas, the array gain allows 
for decreasing the radiated signal energy while maintaining 
the same ASR. However, the energy consumed by circuits 



increases. Maximizing the EE is thus a complicated problem 
where several counteracting factors need to be balanced. This 
stands in contrast to maximization of the ASR, which is rather 
straightforward since the sum capacity is the fundamental 
upper bound. 

There are only a few works in the D2D communication 
literature where the base stations have multiple antennas d)- 
p7) . In | |T3) , uplink MU-MIMO with one D2D pair was 
considered. Cellular user equipments (CUEs) were scheduled 
if they are not in the interference-limited zone of the D2D user. 
The study compared different multi-antenna transmission 
schemes. In |15|, two power control schemes were derived in a 
MIMO network. Two works that are more related to our work 
are and GD- The former investigates the mode selection 
problem in the uplink of a network with potentially many 
antennas at the BS. The impact of the number of antennas on 
the quality-of-service and transmit power was studied while 
users need to decide their mode of operation (i.e., D2D or 
cellular). The latter study, GZ)’ only employs extra antennas 
in the network to protect the CUEs in the uplink. 

The ASR in D2D communications is mostly studied in the 
context of interference and radio resource management GD, 
GD- There are a few works that consider EE in D2D com¬ 
munications, but only for single antenna BSs, e.g., gg, ||2T), 
and | |22) , where the first one proposed a coalition formation 
method, the second one designed a resource allocation scheme, 
and the third one aimed at optimizing the battery life of user 
devices. 

The degrees of freedom offered by having multiple antennas 
at BSs are very useful in the design of future mobile networks, 
because the spatial precoding enables dense multiplexing 
of users with little inter-user interference. In particular, the 
performance for cell edge users, which have almost equal 
SNR to several BSs, can be greatly improved since only the 
useful signal are amplified by precoding ||2^-p5). In order 
to model random user positions, we use mathematical tools 
from stochastic geometry |26|. We assume that there is a fixed 
number of randomly distributed CUEs in the network, while 
the locations of D2D transmitters are randomly distributed 
according to a homogeneous Poisson point process (PPP). We 
characterize the relation between the ASR and EE metrics in 
terms of number of antennas and D2D user density with fixed 
number of CUEs for a given coverage area. 


II. System Model 

We consider a single-cell scenario where the BS is located in 
the center of the cell and its coverage area is a disc of radius R. 
The BS serves [4 single-antenna CUEs which are uniformly 
distributed in the coverage area. These are simultaneously 
served in the downlink using an array of Tc antennas located 
at the BS. It is assumed that 1 < 14 < T). so that the 
beamforming can control interference p7) . 

In addition to the CUEs, there are other single-antenna users 
that bypass the BS and communicate pairwise with each other 
using a D2D communication mode. The locations of the D2D 
transmitters (D2D Tx) are modeled by a homogeneous Poisson 


out-of-cell D2D pair 



Fig. 1. System model where a multi-antenna BS communicates in the 
downlink with multiple CUEs, while multiple user pairs communicate in D2D 
mode. The CUEs are distributed uniformly in the coverage area and the D2D 
users are distributed according to a PPP. The D2D users that are outside the 
coverage area are only considered as interferers. 

point process (PPP) <I> with density in i.e., the average 
number of D2D Tx per unit area is Xd- The D2D receiver 
(D2D Rx) is randomly generated in an isotropic direction with 
a fixed distance away from its corresponding D2D Tx similar 
to the model considered in | |28t . The system model is shown in 
Pig. G] Let Rkj denote the distance between the j-th D2D Tx 
to the fc-th D2D Rx. The performance analysis for D2D users 
is carried out for a typical D2D user, which is denoted by 
the index 0. The typical D2D user is an arbitrarily D2D user 
located in the cell and its corresponding receiver is positioned 
in the origin. The results for a typical user show the statistical 
average performance of the network pS) . Therefore, for any 
performance metric derivation, the D2D users inside the cell 
are considered and the ones outside the cell are only taken 
into account as sources of interference. Note that we neglect 
potential interference from other BSs and leave the multicell 
case for future work. We assume equal power allocation for 
both CUEs and D2D users. Let Pc denote the total transmit 
power of the BS, then the transmit power per CUE is The 
transmit power of the D2D Tx is denoted by Pd- 

Let hj G be the normalized channel response 

between the BS and the j-th CUE, for j G (0,..., Uc — 1}. 
These channels are modeled as Rayleigh fading such that 
hj ^ CA/'(0,1), where CJV{-, •) denotes a circularly symmetric 
complex Gaussian distribution. Linear downlink precoding is 
considered at the BS, based on the zero-forcing (ZP) scheme 
that mitigates interference between the CUEs | [27) . The pre¬ 
coding matrix is denoted by V = [vq, ..., vc/c_i] £ 
in which each column Vj is a normalized precoding vector that 
is assigned to the CUE j in the downlink. Let fo.Bs £ 
be the channel response from the BS to D2D Rx and let 
it be Rayleigh fading as fp^BS ~ CA/^(0,I). Moreover, let 
rj G C and s £ denote the transmitted data signals 

intended for a D2D Rx and the CUEs, respectively. Since each 
user requests a different data, the transmitted signals can be 
modeled as zero-mean and uncorrelated with E[|rjp] = Pd 
and E[||s|p] = Pc- The fading channel response between the 
j-th D2D Tx and the fc-th D2D Rx is denoted by gkj G C 
where gkj ~ CAf{0, 1). Moreover, i?o.BS denotes the random 
distance between the typical D2D Rx and the BS. The pathloss 
is modeled as with i G {c,d}, where index c 

indicates the pathloss between a user and the BS and index d 
corresponds to the pathloss between any two users. Aq^ and 
ai are pathloss coefficient and exponent, respectively. Perfect 






instantaneous channel state information (CSI) is assumed in 
this work for analytic tractability, but imperfect CSI is a 
relevant extension. The received signal at the typical D2D Rx 
is 


yd,0 — \/50,0^'0 + \/^0,cRo BS^ 

'--v-' 

Interference from the BS 

+ \/ ■‘^0,d Ro,i 90,]''’j +'ld,0: 

j^O 


( 1 ) 


Interference from other D2D users 


where r]d,o is zero-mean additive white Gaussian noise with 
power Nq = NqBuj, Nq is the power spectral density of 
the white Gaussian noise, and B^j is the channel bandwidth. 
Therefore, the signal-to-interference-plus-noise ratio (SINR) at 
the typical D2D Rx is 


SINRd,o = 


PdRo,r\9o,o\ 

IbS, 0 + ld,0 + 


No 

Ao,d 


( 2 ) 


in which both numerator and denominator have been normal¬ 
ized by and 


/BS,o = ^^i?-“|||fo%V||2, (3) 

Id.o = Y.PdRoT\9o,]\\ (4) 

3^0 


where /bs,o is the received interference power from the BS 
and Id,o is the received interference power from other D2D 
users that transmit simultaneously. 

Let and eo,fc G C with eo,k ~ CA/'(0,1) be the distance 
and fading channel response between a typical CUE and the 
fc-th D2D Tx, respectively, and let L>o,bs denote the distance 
between a typical CUE and the BS. Then, the received signal 
at the typical CUE is 

9c,0 = 

+ y/AM ^ Dg -fPc.O) (5) 

j 

"-V-" 

Interference from all D2D users 


where pc,o is zero-mean additive white Gaussian noise with 
power Nq. Using the notation C, = Aq^c^, the corresponding 
SINR for the typical CUE is 


SINR,,o = 



Ihp^voP 


No C 
^0,£i ' 


( 6 ) 


where 

= (7) 

3 


is the received interference power from all D2D users. 


III. Performance Analysis 

A. Metrics 

In this paper, two main performance metrics for the net¬ 
work are considered: the average sum rate (ASR) and energy 


efficiency (EE). The ASR is obtained from total rates of both 
D2D and cellular users as 


ASR = UoRc + nR^XdRd, (8) 


where ttR^Xj, is the average number of D2D users in the cell 
and Rt with t G {c,d} denotes the average rates of the CUEs 
and D2D users, respectively. To compute the ASR, we use the 
following lower bound on the average rate for both cellular 
and D2D users |29 Corollary 2] 

Rt = sup B^ log 2 (1 + A) PrjSINRt > /3t}, (9) 

A>o 

which corresponds to the successful transmission rate, and 


= Pr{SINRt > Pt} 


( 10 ) 


is the probability of coverage when the received SINR is 
higher than a specified threshold Pt needed for successful 
reception. Note that SINR^ contains random fading and ran¬ 
dom user locations. Einding the supremum guarantees the best 
rate for the D2D user or CUE. If we know the probability of 
coverage, the expression in (|^ can easily be computed by 
using line search for each user type independently. 

The EE is defined as the ratio between the ASR and the 
total consumed power, i.e.. 


EE = 


ASR 

Total power 


( 11 ) 


Eor the total power consumption, we consider a detailed model 
described in j?) as 

Total power = - (Pc + Xd^R^Pd) + Cq 

Tj ' 

+ TeCi + (Uc + 2Xd^R^)C2, (12) 


where Pc + Xd^R^Pd is the total transmission power, rj is the 
amplifier efficiency (0 < ry < 1), Cq is the load independent 
power consumption at the BS, Ci is the power consumption 
per BS antennas, C 2 is the power consumption per handset, 
and Uc + 2XdTR^ is the number of active users. 


B. Probability of Coverage 

In order to calculate the ASR and EE, we derive the 
probability of coverage for both cellular and D2D users. This 
is one of the main contributions of the paper. Due to the page 
limit, we only provide outlines of the proofs. 


Proposition 1. The probability of coverage for a typical D2D 
user is given by 


-pd2d _ 
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0,d 


, (13) 





















where x = 

sinc(z) = 
function. 
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and B{y, a, h) is the incomplete Beta 


Proof: This result follows by substituting the definition 
of SINRd 0 from (|^ into ( [TOl i where we obtain 

Rni 


pd2d _ p 
^succ ^^ 


or > 


m,o 


W IP 

- ^7BS,o.7d.( 


Pd 

— Pdd:{lBS,0+td,0 + f^ 


/3d(^BS,0 + ld,0 + I 


- PlBS,oi^dx)Ci^^^{Pdx)e 




(14) 


where C denotes the Laplace transform defined as Cz{s) = 
Step (a) comes from the fact that |po,or ’’’’ exp(l) 
and (6) follows since the interference terms are independent 
and the noise term is not dependent on the interference terms. 
The first Laplace transform in ( [l4| is with respect to /bs.o 
in which is a function of two random variables, i.e., 
Ilfr^gjVir where the norm is well-approximated by a Chi- 


squared distribution as 2 ||fr^g 5 V|p ^ xic/c ^o.bs 


which is uniformly distributed over the coverage area. Cal¬ 
culating this Laplace transform results in the non-exponential 
term of ( [T3| ). The second Laplace transform in ( [T4l l is with 
respect to Id,o in •Hji and is calculated based on the probability 
generating functional (PGFL) |26| which results in the first 
exponential term in ([T3|l. ■ 


Proposition 2. The probability of coverage for a typical 
cellular user is given by 


Kucc = Edo., 


^ An J 


Tc-t/c 


k—0 i—0 




(-1) TrexP —^ 


ds* 


sinc( —) 


(15) 


where s = ^A^^dDoBs ^ stands for the i'^ derivative. 
Proof: Substituting SINRcp into ( [T0| ), we get 

PLc = Pr jlh^vop > ^Ao,dD^:^,[ld,a + 


(a 


E 


Do^BsJd^i 




Tc-C/c 




Alo,d 

No \>^ 


k=0 


A, 


0,d 


(16) 


where (a) follows from the CCDF of jh^vop with 
2|h^VoP - X 2 (rc-( 7 c+i) ^o,BS and Id,c- Taking the 

expectation with regard to the interference Id^c gives m 

The expectation in ( [T5] l with respect to Hq.bs can be com¬ 
puted numerically. The analytical results of Proposition [T] and 
Proposition have been verified by Monte-Carlo simulations. 
A main benefit of the analytic expressions is that they can 
be computed efficiently, which basically is a prerequisite for 
the multi-variable system analysis carried out in Section 
Using the results from Proposition [T] and Proposition we 


TABLE I 

System and simulation parameters. 


Description 

Parameter 

Value 

D2D TX power 

Pd 

13 dBm 

BS TX power 

Pa 

41 dBm 

Number of CUE 

Ua 

4 

Cell radius 

R 

500 m 

Bandwidth 


20 MHz 

Thermal noise power 

No 

-131 dBm 

Noise figure in UE 

F 

5 dB 

Carrier frequency 

fa 

2 GHz 

D2D pair distance 

ddid 

50 m 

Pathloss exponent betw. devices 


3 

Pathloss exponent betw. BS-device 

CAc 

3.67 

Pathloss coefficient betw. devices 

Ad 

38.84 dB 

Pathloss coefficient betw. BS-device 

Aa 

30.55 dB 

Amplifier efficiency 

V 

0.3 

Load-independent power in BS 

Co 

5 W 

Power per BS antenna 

Cl 

0.5 W 

Power per UE handset 

C2 

0.1 W 


proceed to evaluate the network performance in terms of the 
ASR and EE from (|^ and ( [TT] l, respectively. 

IV. Numerical Results 

In this section, we assess the performance of the setup in 
Fig-0 in terms of ASR and EE using numerical evaluations. 
Both EE and ASR are functions of three key parameters, 
namely, the number of BS antennas Tc, the density of D2D 
users Xd, and the number of cellular users Uc. Due to space 
limitations, we cannot show the individual effect of all param¬ 
eters. Therefore, we fix the number of CUEs to (7c = 4 and 
study the tradeoffs among other parameters. The system and 
simulation parameters are given in Table 

In Eig. 1^ the behavior of ASR is shown with respect to 
different number of antennas Tf. and the density of D2D users 
Xd. It is observed that increasing the number of antennas 
always seems to increase the ASR. In contrast, for the ASR, 
there is an optimal value of the D2D density Ad, which results 
in the maximum ASR for each number of antennas. However, 
there is a difference in the shape of the ASR in lower Tc values 
and higher Tc values. In order to clarify this effect, we plot 
the ASR versus Xd in a 2-D plot with Tc = {4, 70} in Eig. 0 

Eor Tc = 4, the rate contributed from CUEs to the sum rate 
is low and adding D2D users to the network (i.e. increasing 
Xd), even though they create interference, leads to increasing 
ASR until a certain density. This increase in the ASR continues 
until a point where the cross-tier interference between D2D 
users limits per link data rate and results in decreasing ASR. 
By increasing the number of antennas to Tc = 70, the rates 
of the cellular users becomes higher. However, by introducing 
small number of D2D users, the effect of the initial interfer¬ 
ence from D2D users becomes visible, i.e. the decrease in the 
CUEs’ rates is not compensated in the ASR by the contribution 
of the D2D users’ rate. However, if we keep increasing the 
number of D2D users, even though the rate per link decreases 
for both CUE and D2D users, the resulting aggregate D2D 
rate becomes higher and the ASR starts to increase, which 
is the reason behind the local minima. The second turning 
point follows from the same reasoning as with Tc = 4, i.e. in 
higher D2D densities, the interference from D2D users are the 
limiting factor for the ASR. This effect can also be observed in 
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Fig. 2. ASR [Mbit/s] as a function of D2D user density and BS antennas 

Tc. 

Fig. 1^ where the ASR performance is depicted versus different 
number of antennas for two D2D densities. At lower densities, 
increasing the number of antennas is beneficial in terms of 
ASR, however, in the interference-limited regime (higher Ad), 
increasing the number of antennas does not impact the total 
network performance. 

Fig.|^shows the network performance in terms of EE versus 
the parameters Xd and Tc- To study this results further, similar 
to the ASR, we first plot the EE versus Ad for Tc = {4, 70} 
in Eig. 1^ We can see that the pattern for both lower and 
higher number of antenna is similar to Fig. However, if we 
plot the EE versus Tc, we see a different behavior for low- 
and high-density D2D scenario. Erom Eig. it is observed 
that indeed the low-density scenario benefits in terms of the 
number of antennas until the sum of the circuit powers for all 
antennas dominates the performance and leads to a gradual 
decrease in EE. As the figure implies, there exists an optimal 
number of antennas (Tc) for maximal EE in the low-density 
scenario. However, in high-density D2D scenario, which is an 
interference-limited scenario, the EE decreases almost linearly 
with Tc, as depicted in Eig. Increasing the number of 
antennas in this region cannot improve the ASR much as 
shown in Fig. Then, at the same time the circuit power 
resulting from higher number of antennas increases which in 
turn leads to poor network EE. 

V. Conclusions 

We studied the coexistence of two key 5G concepts: device- 
to-device (D2D) communication and massive MIMO. We 
considered two performance metrics, namely, average sum rate 
(ASR) in bit/s and energy efficiency (EE) in bit/Joule. We 
assumed that there is a fixed number of randomly distributed 
cellular user equipments (CUEs) in the network, while the 
number of D2D transmitters are distributed according to a 
Poisson point process (PPP). We derived tractable expressions 
for ASR and EE, and studied the tradeoff between the number 
of base station antennas and the density of D2D users with 
a fixed number of CUE for a given coverage area in the 
downlink. Our results showed that both ASR and EE have 



Fig. 3. ASR [Mbit/s] as a function of D2D user density for Tc = 
{4,70}. 



Fig. 4. ASR [Mbit/s] as a function of BS antennas Tc S (4,..., 100} for 
Ad = {10-6,10-4}. 

different behaviors in scenarios with low and high density 
of D2D users. By increasing the number of antennas in the 
low D2D density regime, the ASR improves, however, EE 
increases until the circuit power from many antennas becomes 
dominant. In the high D2D density regime, there is small 
gain in terms of the ASR from adding many antennas to the 
detriment of EE which decreases tremendously. 
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